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In this paper, we employed a non-critical QCD holographic model constructed in the AdS& su- 
pergravity to study a baryon. We show that the size of baryon is of order one respect to the A but 
also is smaller than the scale of the dual QCD. We obtain an effective action for the four dimen- 
sional nucleon in terms of the meson exchange potentials and calculate all meson-nucleon couplings 
in AdSa background. We consider four modern phenomenological interaction models and compare 
our results with their predictions. Also, our numeric results are compared with the results of the 
Sakai-Sugimoto model results which indicate that the non-critical holographic QCD model can be 
a good toy to calculate the meson-nucleon couplings. 
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I. INTRODUCTION 

Construction of nucleon-nucleon (NN) potential has a 
long history in nuclear physics due to its role in under- 
standing the nuclear force. Many potential models have 
been constructed from the 1950s which have been sup- 
posed to fit the available NN scattering data. The newer 
potentials have only slightly improved respect to the pre- 
vious ones in describing the newer and much more accu- 
rate data. As it is shown in ref[l], all of these potential 
models don't have a good quality with respect to the pp 
scattering data below the 350 MeV and just a few of them 
are of a satisfactory quality. These models are the Reid 
soft-core potential [2] Reid68, the Nijmegen soft-core po- 
tential [3] Nijm78, new Bonn pp potential [4] Bonn89 and 
also the parameterized Paris potential [5] Paris80. These 
familiar one-boson-exchange potentials (OBEP) contain 
a relatively small number of free parameters (about 10 to 
15 parameters), but haven't a reasonable description of 
the empirical scattering data. Also, most of them which 
have been fitted to the np scattering data unfortunately 
do not automatically fit the pp data even by considering 
the correction term for the Coulomb interaction [1]. Of 
course, new versions of these potentials have been con- 
structed such as Nijm I, Nijm II, Reid93 [6], CD-Bonn [7] 
and AV18 [8] which explain the empirical scattering data 
successfully. But they contain a large number of purely 
phenomenological parameters. For example, an update 
(Nijm92pp [9]) of the Nijm78 potential contains 39 free 
parameters which determined using the scattering data. 

On the other hand, there are many attempts to impose 
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the symmetries of QCD using an effective Lagrangian of 
pions and nucleons [10,11]. These models only capture 
the qualitative features of the nuclear interactions and 
couldn't compete with the much more successful poten- 
tial models mentioned above. Despite many efforts, no 
potential model has not yet been constructed which give 
a high-quality description of the empirical data, obey the 
symmetries of QCD and contain only a few number of free 
phenomenological parameters. 

One of the applications of AdS/CFT duality [12-14] 
is holographic QCD introduced recently to solve the 
strong-coupling QCD problems such as chiral dynam- 
ics of hadrons in particular baryons [15-38]. Among the 
holographic QCD models, it can be noted to the Sakai- 
Sugimoto (SS) [39, 40] and Klebanov-Strassler (KS) mod- 
els [41]. 

The predictions of the SS model are in a good agree- 
ment with the lattice simulations such as a glueball spec- 
trum of pure QCD [42, 43]. Also this model describes 
baryons and their interactions with mesons [22-24, 39, 
40]. It is shown that the baryons can be taken as the 
point-like objects at distances larger than their sizes, so 
their interactions can be described by the exchange of 
light particles such as mesons. Therefore, one can find 
the baryon-baryon potential from the Feynman diagrams 
by using the interaction vertices including baryon cur- 
rents and light mesons [23]. But there are some incon- 
sistencies. For example, the size of the baryon is pro- 
portional to A -1 / 2 . Consequently in the large 't Hooft 
coupling (large A) the size of the baryon becomes zero 
and the stringy corrections have to be taken into account. 
Another problem is that the scale of the system associ- 
ated with the baryonic structure is roughly half the one 
needed to fit to the mesonic data [44]. 

Also there is another problem for such holographic 
models arising from the critical string theory. In these 
models, the color brane backgrounds are ten-dimensional 
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so the dual gauge theories are super symmetric. In or- 
der to break the super symmetry, some parts of such 
backgrounds need to be compactified on some manifolds. 
This causes to produce some Kaluza-Klein modes with 
the mass scale of the same order as the masses of the 
hadronic modes. These unwanted modes are coupled to 
the hadronic modes and there is no mechanism to dis- 
entangle them from the hadronic modes yet. In order 
to overcome this problem, it is possible to consider the 
color branc configuration in the noncritical string theory. 
The result is a gravitational backgrounds located at the 
low dimensions [45-48]. In these backgrounds the string 
coupling constant is proportional to so the large N c 
limit corresponds to the small string coupling constant. 
However, contrary to the critical holographic models, in 
large N c limit the 't Hooft coupling is of order one in- 
stead of infinity and the scalar curvature of the gravita- 
tional background is also of order one. So, it seems the 
noncritical gauge-gravity correspondence is not very re- 
liable. But studies show that the results of these models 
for some low energy QCD properties such as the meson 
mass spectrum, Wilson loop and the mass spectrum of 
glueballs [49-51] are comparable with the lattice compu- 
tations. Therefore, the noncritical holographic models 
seem still useful to study the QCD stuff. 



II. HOLOGRAPHIC QCD FROM 
NON-CRITICAL STRING THEORY 

In the presented noncritical model, the gravity back- 
ground is generated by near-extremal DA branes wrapped 
over a circle with the anti-periodic boundary conditions. 
Two stacks of flavor branes, DA branes and anti-D4 
branes, are added to this geometry which called flavor 
probe branes. The color branes extend along the direc- 
tions t, Xi,X2,Xs,t while the probe flavor branes fill the 
whole Minkowski space and stretch along the radius U 
up to infinity. The strings attached color D4-brane to a 
flavor brane transform as quarks, while strings hanging 
between a color DA and a flavor DA transform as anti- 
quarks. The chiral symmetry breaking is achieved by a 
reconnection of the brane-anti-brane pairs. Under the 
quenched approximation (N c ^> Nf), the backreactions 
of flavor branes on the color branes can be neglected. 
Just like the SS model, the t coordinate is wrapped on 
a circle and the anti-periodic condition is considered for 
the fermions on the thermal circle. The final low energy 
effective theory on the background is a four-dimensional 
QCD-like effective theory with the global chiral symme- 
try U(N f ) L x U(N f ) R . ' 

In this model, the near horizon gravity background at 
low energy is [49] 



One of the noncritical holographic models is composed 
of DA and anti DA brane in six-dimensional noncritical 
string theory [47,49]. The low energy effective theory on 
the intersecting brane configuration is a four-dimensional 
QCD-like effective theory with the global chiral symme- 
try U(Nf)L x U(Nf)n. In this brane configuration, the 
six-dimensional gravity background is the near horizon 
geometry of the color DA brane. This model is based on 
the compactified AdS§ space-time with a constant dila- 
ton. So the model does not suffer from the large string 
coupling as the SS model. The meson spectrum [49] and 
the structure of thermal phase [52] are studied in this 
model. Some properties, like the dependence of the me- 
son masses on the stringy mass of the quarks and the 
excitation number are different from the critical holo- 
graphic models such as the SS model. 
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where f(U) and RR six- form field strength F( 6 ) are de- 
fined by the following relations, 
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For the solutions of near extremal flavored AdSe, the 
values of dilaton and Ra<is are as follows, 



This paper is organized as follows: In section 2 we 
briefly review the noncritical model and the mode ex- 
pansion of the gauge field there. In section 3 we analyze 
the baryon and extract its mass and size. We derive the 
effective action for the baryon and obtain a noncritical 
prescription of nucleon-nucleon potential in terms of the 
meson exchange interactions in section 4. In section 5 
we obtain the nucleon-meson couplings and compare our 
results with the some available data predicted in four 
modern phenomenological models (Nijmegen (93), Paris, 
CD-Bonn and AV 18 models). Section 6 is devoted to 
brief summary and conclusions. 
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This relation indicated that the Ra<is and dilaton depend 
on the ratio of number of colors(^ Q c ) and flavors(~ Qf). 
Under the quenched approximation that we considered 
here, the values of the dilaton and AdS radius become, 



R 2 



AdS 



15 

y 



2V2 

V3Q C 



(4) 



3 



where Q c is proportional to the number of color branes, 
N c . 

To avoid the singularity, the coordinate t satisfies the 
following periodic condition, 

AttR 2 



t <~ r + St , 



St 



(5) 



5U KK 

Also, the Kaluza-Klcin mass scale of this compact dimen- 



sion is 
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M KK = — = 



KK 



St - 2 # ' (6) 

and dual gauge field theory for this background is non 
supersymmetric. 

In AdS/QCD we have gauge field living in the bulk 
AdS whose dynamics is dual to the meson sector of QCD, 
including pions and higher resonances. The gauge field 
on the DA brane has five components, A ll (fi = 0, 1, 2, 3) 
and Ajj. The DA brane action is given by 

Sda = -M4 J d 5 xpe~ 4> y / - dct(g M N + 2-Ka'F M N) 

+S C s, (7) 

where the second term is the Chern-Simons action and 
yU4 = 2tt/(2ttIs) 5 . It is useful to define the new variable 
z as follows 



u, = m K + u i KK z i ) 



2x1/5 
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Then by neglecting the higher order of F 2 , the DA brane 
action becomes 

R 4 
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where we define 

« = 4^f^ (10 » 

The gauge fields (/i — 0,1,2,3) and A z can be ex- 
panded in terms of complete sets {ip n (z)} and {<fi n (z)} 
as 

A^,z) = Y,Bf\x*)M*) > (11) 

n 

A z {x»,z) = 5><"V)0n(*) ■ (12) 

n 

By computing the field strengths, the action (9) is written 

as 

R 4 



'DA 



-^ 4 (2W) 2 J d 4 xdz 



Let us consider just the vector meson field B { ™ ] . So, we 
need to keep the following part of action 
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We introduce the dimcnsionlcss parameters as follows 
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(14) 



K(z) = l + z = 



~2 ( U z 
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Ukk \Ukk 
and rewrite the action (14) in terms of these quantities 
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We choose the normalization condition for ip n (n > 1) as 
R 4 
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Also, we suppose the functions tp n (n > 1) satisfy the 
following condition 

M 4 (2W) 2 4t 2 - / dzK 9 / 10 9^ m &^ n = X n 8 nm (18) 
Ukk J 

We find an eigenvalue equation for the functions ip n (n > 
1) from the equations (17) and (18) as 



-K^ 2 & ; (K^ c^p m )=X m ^ 



(19) 



Considering the above conditions, action becomes canon- 
ically normalized 



Sd4 
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where is a massive vector meson of mass m„ = 

\]J 2 Mkk for all n > 1. Let us include tpW and rewrite 
the pseudoscalar part of action (13) in terms of new vari- 
ables (15) 

Sda = -M 4 (2W) 2 j d 4 xdI^-U 3 K / K K g / 10 J2 
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p In order to normalize the kinetic part of the above 
action, we consider the ortho-normal condition for ip( n ) 
as 



,<t>n) = — M 27ra ) U^k 
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(22) 



If we multiply the equation (22) by A„ and compare 
it with equation (18), it is obvious that the functions 
and ip n relate together. In fact, we can consider 
4> n = m^ 1 ip n (n > 1). Also, there exist a function 
4> = C/K 9 / 10 that is orthogonal to ip n for all n > 1 



{4>Q,4> n ) oc J dzd^ip n = 0. (for n > 1) (23) 

We can use the normalization condition 1 = (</>o,</>o) to 
obtain the normalization constant C. By using an ap- 
propriate gauge transformation, the action (9) becomes 



~>D4 



d 4 x 
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is the pion held, which is the Nambu-Goldstone bo- 
son associated with the chiral symmetry breaking. In- 
terpretation of this held is like the one in critical models 
such as the SS model [39]. Therefor it is not necessary 
to repeat it here. 

Usually it is useful to make another gauge choice, the 
A z = gauge, to be sure that the field strengths are 
vanish at z — > ±oo. Actually we can transform to the 
new gauge through the gauge transformation 

A M -> A M - d M A , (25) 

and obtain the following new gauge fields 

A z (x^z) = , 

A^,z) = -9^°V)V>o(z) + E4"V)V>n(2) • 

n>l 

(26) 

Function ipo(z) is calculated as 



= / dz'Mz') = CU KK zF 1 (0.5,0.9,1.5,-z 2 ). 
Jo 

(27) 

where F\ is a Hypergeometric function. It is noted that 
the massless pseudo-scalar meson appears in the asymp- 
totic behavior of A^ since we have 

A^.z) -> ±1.8CU KK d^°\x») (as z -> ±oo).(28) 

In order to obtain the meson spectrum, it is needed to 
solve the equation (19) numerically by considering the 
normalization condition (17). 



Now, we have to solve (19) with the normalization 
condition given by (17) numerically to obtain the meson 
spectrum. Since Eq. (19) is invariant under z — > —z we 
can assume ip n to be an even or odd function. In fact, the 
B^ is a four-dimensional vector and axial vector when 
ip n is an even or odd function, respectively. We solve the 
Eq. (19) by means of shooting method and obtain the 
mass of lightest mesons. We compare our results with 
the SS model results and experimental data in Table. 1. 
As you can see, our result are in a good agreement with 
the previous ones. Also, the same results obtained in the 
[49] in the AdSe background which is exactly coincident 
with our results. 



TABLE I: . The ratio of the obtained eigenvalue compared 
with the SS model results[39] and the ratio of meson masses. 
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Experiment 
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1.8 


m 2 
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We can generalized above analysis to the case of Nf > 
1 flavor QCD by introducing Nf probe D4-branes. In 
order to obtain a finite four-dimensional action for the 
modes localized around z = 0, the field strength Fmn 
should vanish at z = ±oo. This implies that the gauge 
field Am must asymptotically take a pure gauge config- 
uration 



A M (x M , z) -> U± 1 (x ti , z)d M U±(x' 1 , z) . (as z -> ±oo) 



(29) 



In analogy with the analysis done in SS model [39], we 
can write 



A^ X ^z) = U-\x^U(x^ + (z) + Y J B { ;\x^ n (z) . 

n>l 

(30) 

where 



(31) 



^o(2)=3g ^i(0.5,0.9, 1.5,- F). (32) 

Now, let us neglect the vector meson fields B^ ] (n > 1) 
and write the field strengths as 



F^ = [U^d^U-^U] ^+(^+-1) 
F ZI1 = U^dpU&jioiz), 



(33) 
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Substituting these quantities into the non-Abelian gen- 
eralization of action (9), we obtain 

S D 4= - M^a' f j d 4 x tr (^{U^d^Uf 

+ B[U-%U,U- 1 d v U] 2> j , (34) 

where the coefficient A and B arc defined as follows 



ok tj 1 / 2 

2 _ zo u KK 

A 3.6 ' 



B = *Tj dZ W^-V 2 = W^- (35) 



If we compare the action (34) with the familiar Skyrme 
model action 



S = J d^x^tviU^d^U) 2 



(36) 



it is possible to compute the pion decay constant and 
dimcnsionless parameter e in terms of the noncritical 
model parameters 



/^ = 4m 4 (2W)^=W- 



3 45/^(2™')' 
3.6 R 3 



N C M 2 KK (37) 



and 



\ = 32^ 4 (2W) 2 S = J Indira') 2 RN C (38) 



It is clear from the above equations that parameters /„■ 
and e depend on N c as ~ 0(y/NZ) and e <~ 0{l/\fWl) 
respectively, which are in agreement with the result ob- 
tained from the SS model and also QCD in large N c . We 
fix the M KK such that the f„ ~ 93 McV for N c =3. So, 
we obtain Mkk — 395 MeV for our holographic model. 
It should be noted that Mkk is the only mass scale of 
the model. 



III. BARYON IN AdS 6 

In this section we aim to introduce the baryon config- 
uration in the non-critical holographic model. As it is 
known, in the SS model the baryon vertex is a DA brane 
wrapped on a S 4 cycle. Here in six dimension, there is no 
compact S* 4 sphere, so we introduce an unwrapped DO 
brane as a baryon vertex instead. In analogy with the 
SS model, there is a Chern-Simons term on the vertex 
worldvolumc as 



Scs oc 



J dtA (t), 



(39) 



which induced N c units of electric charge on the un- 
wrapped DO brane. In accordance with the Gauss con- 
straint, the net charge should be zero. So, one needs 
to attach N c fundamental strings to the DO brane. In 
turn, the other side of the strings should end up on the 
probe DA branes. The baryon vertex looks like an ob- 
ject with N c electric charge respect to the gauge field on 
the DA brane whose charge is the baryon number. This 
DO brane dissolves into the DA brane and becomes an 
instanton soliton. It is important to know that what is 
the size of instanton in our model. In the SS model it is 
shown that the size of such instanton goes to zero at large 
't Hooft coupling constant which is one of the problems 
of the SS model in describing the baryons [23]. 

Let us consider the DBI action in the Yang-Mills (YM) 
approximation for the DA brane 



>YM 



i^ 4 (2W) 2 J e ^ v /q^i"trF mn F mn . (40) 



The induced metric on the DA brane is 
2 / 



54+1 



Now, we define the new coordinate w as 

R 2 U l ' 2 dU 

dw = 



(42) 



Using this new coordinate, the metric (41) transforms to 
a conformally flat metric 



where 



54 + i = H(w) [dw 2 + rj^dx^dx") 



H(w) = (|) 2 . 



(43) 



(44) 



Also, we can rewrite the w coordinate in terms of equa- 
tion (15) as 



2 R 2 Ul K dz 
dw = - - 



(45) 



Note that in the new conformally flat metric, the fifth 
direction is a finite interval [—w ma x,w m ax\ because 



00 R 2 U^ 2 dU 
R 2 r 00 dU R 2 



r z r 

UkkJi 



Ukk Jl ^JJ5 - 1 Ukk 
We can approximate w near the origin w ~ 0, as 



2 / R 
w ~ - — — I z. 



1.25 < oo.(46) 



5 \ Ukk 



(47) 
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and using the relation (6), it is rewritten as 



The size of an stable instanton comes from the minimiz- 
ing the total energy 



w 



M KK U KK 
or equivalently, 



or Mkk w ~ 



Ukk 



U^U 5 KK (l + M] (K w 2 ) 



(48) 



(49) 
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3/2 27^(2^ a') 2 M 2 KK ' 
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This indicates that, here in analogy with the SS model, 
Mkk is the only mass scale that dictated the deviation 
of the metric from the flat one and it is the only mass 
scale of the theory in the low energy limit. (Please note 
that the DA branes come with two asymptotic region 
at w — > ±w max corresponding to the UV physics and 
infrared region near the w ~ 0. ) 

The action (40) is rewritten in conformally flat metric 
(43) as 



-.154 



~ P4(2ira') 2 J d 4 xdwe 



U{w) 
R 



trF mn F mn 



/ 



= / dx 4 dw 



trF F mn 

Ae z (w) 



(50) 



So, the position dependent electric coupling e(w) of this 
five dimensional YM equals to 



1 _ ^.^aVRN, ^ MM 
e z (w) 5 \UkkJ 

As it is known for a unit instanton 

J trF A F = 1 J trF mn F mn = 8tt 2 , (52) 

Inserting the above relation in the equation (50), we 
obtain the energy of a point-like instanton localized at 

w = as 



In the SS model (the critical version of dual QCD) the 
size of instanton goes to zero because of the large 't Hooft 
coupling limit. But in noncritical string theory, the 't 
Hooft coupling is of order one implying the size of in- 
stanton also be of order 1. But here it is still smaller 
than the effective length of the fifth direction <~ 1 /Mkk 
of the dual QCD. 



IV. NUCLEON-NUCLEON POTENTIAL 

In the previous section, we demonstrated that the size 
of baryon in the non critical holographic model dual to 
QCD is smaller than the scale of the dual QCD and we 
can assume that the baryon is a point-like object in five 
dimensions. Thus as a leading approximation, we can 
treat it as a point-like quantum field in five dimension. 
In this paper, we will restrict ourself to the fermionic 
baryon because we intend to study the nucleons. So, we 
must consider the odd N c to have a fermionic spin 1/2 
baryon. Also, we will assume Np = 2 and consider the 
lowest baryons which form the proton-neutron doublet 
under SU(Np = 2). All of these lead us to introduce 
an isospin 1/2 Diprac field M for the five dimensional 
baryon. 

The leading 5D kinetic term for N is the standard 
Dirac action in the curved background along with a po- 
sition dependent mass term for the baryon. Moreover, 
there is a coupling between the baryon field and the gauge 
filed living on the flavor branes that should be considered. 
Therefore a complete action for the baryon reads as 



m 



(0) 



3/2 47H//4 (2na'YR 



N c M KK . (53) 



d xdw 



iAfj m D m Af - im b {w)MN 



By increasing the size of instanton, more energy is needed 
because l/e 2 (i«) is an increasing function of \w\. So in- 
stanton tends to collapse to a point-like object. On the 
other hand, N c fundamental strings attached to the DA 
branes behave as N c units of electric charge on the brane. 
The Coulomb repulsion among them prefer a finite size 
for the instanton. Therefor, there is a competition be- 
tween the mass of instanton vs. Coulomb energy of fun- 
damental strings. For a small instanton of size p with the 
density D(x l ,w) ~ p 4 / (r 2 + w 2 + p 2 ) 4 , the YM energy is 
approximated 



+ 35M 



Pba 



ryon R~i mn F mn M 



e 2 (w) 



1 



— Id xdw 



Ae 2 (w) 



trF mn F D 



(57) 
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^-m {0) M 2 n 2 
g m B 1V1 KKP ' 

and the five dimensional Coulomb energy goes as 
1 e(0) 2 A 2 



IOttV 



(54) 



(55) 



where D m is covariant derivative, pbaryon is the size of the 
stabilized instanton and 175 (w) is an unknown function 

that its value at w = is 2tt 2 /3 [23]. The factor p ^°" 
is putted just for some convenience. Usually, the first 
two terms in the action are called the minimal coupling 
and the last term in first integral refers to the magnetic 
coupling. 

Four dimensional nucleons are the localized modes at 
w ~ which are the lowest eigenmodes of 5D baryon 
along the w direction. So, the 5D baryon action must 
reduces to the 4D. We can do it using performance the 
KK mode expansion for the action. In order to do that, 
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we must mode expand the baryon field M{x lll w) and the 
gauge field A{x )1 ,w). The gauge field has a KK mode 
expansion which studied in section (3) in detail. The 
baryon field can be expanded as 



Af L , R (^,w) = N L , R (x»)f L , R (w), 



(58) 



where Nl^ r (x^) is the chiral component of the 4D nu- 
cleon field. The profile functions fL,it(u>) satisfy the con- 
ditions 

d w fL(w) +m b (w)f L (w) = m B jFi{w) , 
-d w f R (w) +mb(w)f R (w) = m B fL{w) , (59) 

in the range w <G [—w max ,w ma x}, and the eigenvalue tub 
is the mass of the nucleon mode N(x). Also, the eigen- 
functions /l,_r(w) obey the following normalization con- 
dition as 



/" 

J — ', 



dw \Jl{w)\ 2 = 



dw \f R (w)\ 2 = l, (60) 



It is more useful to consider a second-order equation for 

Il,r( w ) 

[-dl - d w m b (w) + (m b (w)) 2 ] f L (w) = m%f L (w) , 
[-dl + d w m b (w) + (m b (w)) 2 ] f R (w) = m%f R (w) . 

(61) 

As we approach w — > ±w max , m b (w) diverges as 
m b (w) ~ — p and the above equations have nor- 
malizable eigenfunctions with a discrete spectrum of tub- 
Note that the term —d w m b (w) is asymmetric under 
w — > —w. It causes to /l(w) tends to shift to the posi- 
tive w side, and the opposite happens for f R (w) and it is 
important in axial coupling of the nucleon to the pions. 

As it is indicated in the section (2), the gauge field has 
a mode expansion (30) in the A z = gauge. We can 
rewrite the equation (30) as 

An(x,w) = ia^(x)^o(w) + iP^x) Bl n) (x)ip {n) (w) . 

n 

(62) 

Here, we use the above expansion along with the prop- 
erties of functions fh{w) — ±/r(— w) and V>o and ip(n) 
under w — >• — w to calculate the 4D action. It is worth- 
while to note that again the ip(2k+i){ w ) is even, while 
4 l (2k){ w ) is °dd under w — > —w, corresponding to vector 
B^ k+1 \x tl ) and axial- vector mesons B^ k \x^), respec- 
tively. For simplicity, we neglect the Chern-Simons term 
in the baryon action (57). By inserting the mode ex- 
pansion of the baryon field in the action, we obtain the 
minimal coupling as 

^minimal 

J d A xdw [-iN fj m (d m - iA m )Nf - im b {w)N fNf] 



= J d 4 x [-iN Yd^N - im B NN] 
- J d 4 xdw [Nf^A^Nf] 



Now, we expand the gauge field presented in the last inte- 
gral using the equation (62). Since the parity of ipr n \(w) 
depends on n, we able to separate the odd and even n. 
After taking the integrals over to, we obtain the 4D min- 
imal action for the nucleon as 



(63) 



^minima, 



'=J d ' 



iN'fd^N - im B NN 



nminimaL _■_ nmimmal 
-^vector ' ^axial 



(64) 



where the minimal vector and axial interactions are 

k>0 

(65) 



k>l 



The various minimal couplings constants 9v^ min , g A ] min 
as well as the pion-nucleon axial coupling gA,min are cal- 
culated by suitable wave-function overlap integrals as 



(fe) 

rntn 



dw \f L {w)\ V(2fc+1) (w) 



9 A ,min = I dw \hH\ 1>(2k)(v>) 



9A,r. 



dw |/l(w)| ipo(w) . (66) 



Also, the magnetic interaction term in the action (57) 
becomes 



^magnetic 



lc = - J d 4 xdw 



(67) 



Inserting the gauge field expansion into the equation (67), 
the magnetic interaction reads as 

^magnetic _ J ( jA x ^magnetic _j_ £magnetic\^ (Q8) 



where 



^magnetic 
-'vector 



= -T.^lagNl^B^N 



k>Q 



^magnetic _ WA,mag J\f^^ a ^]\T 



- E^,U^7 5 4 2fc) iV, (69) 



fc>i 



and the magnetic coupling are defined as 



9A 



i nag 



/v 



dw ( ff 5 (w) Pb °7T I l/i( w )| 2 d wi>o{w) 
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(fe) 

,mag 



(fe) 
@V,mag 



/W ma ,x / 
dw g 5 (w) 
-w max \ 



baryon 

e 2 (w) 



|/l(w)| d w ip(2k)(w) 



/V 



dw g 5 (w) 



Pbaryon 

e 2 (w) 



\Jl(w)\ d w tp {2k+1) (w). (70) 

Using equation (51), we can rewrite the magnetic cou- 
plings as 



J*0 on 

9v,mag ~ * ^mag 



n (k) - on 

9A,mag ~ ^^mag 



9A,mag 4 C m ag 



Where we define C mag as 

l\2 



[ Wma ° . (g 5 (w)\ (u{w) 



\!l{w)\ d w 1p(2k+l){w) , 



35(0)7 \ Ukk, 
\Jl{w)\ 2 d w ip(2k)M , 

\ 95(0) J \ Ukk) 



3/2^4 (2 it a') 



C mag = v-i-->y-"-> RNc 6 ( O ) M xx p^on -(72) 

Also, there is a next leading order term in the magnetic 
coupling action which causes the derivative couplings, fi- 
nally by consideration the derivative terms, the nucleon 
Lagrangian is written as 

£ N ucicon = -iN-f^d^N - im B NN + £ vcc tor + £ axial , 

(73) 

and 

fe>0 



fe>0 



fc>l 



(74) 



fe>0 



where g = g m in + g-mag for all couplings. We neglect 
the derivative couplings in the following calculation as 
leading approximation. 

It should be noted that the iso-scalar mesons couple to 
the nucleon in a different formalism than the iso-vector 
mesons, since instanton carries only the non-Abelian held 
strength. So for the iso-scalar mesons, such as lo^ me- 
son, only the minimal couplings contribute 

singlet 

9 A g A .min > 



(k)singlet (fe) 

9 A ~ 9A,min ' 

(k)singlet (fe) 

9v = 9v,min ■ 



(75) 



But the iso-vector mesons couple to the nucleon from 
both minimal and magnetic channels. Thus, the iso- 
vector meson couplings are 



triplet 
9 A 


QA,min 


\~ QA,mag > 




^(k)triplet 


(fe) 

9 A,min 


(fe) 

9 A,mag ' 




(k)triplet 

9v 


9v,min 


(fe) 
9v,mag • 


(76) 



The iso-scalar and iso-vector mesons have the same origin 
in the five dimensional dynamics of gauge field. In fact, 
if we write the gauge field in the fundamental represen- 
tation, we could decompose the massive vector mesons 
as 



B (2fc+i) = f 1/2 \ (fe) (fe) 
^ 1/2 ) » P » 



(77) 



where lo^ and p^ arc the iso-scalar and the iso-vector 
parts of a vector meson. Since baryon is made out of 
N c product quark doublets, the above composition for 
nucleon should be written as 



„(2fe+l) 



N c /2 
N c /2 



u,W+pW. (78) 



Therefor, there is an overall factor N c between the iso- 
scalar ui^ and iso-vector p^ mesons. Indeed, there is a 
universal relation between the Yukawa couplings involv- 
ing the iso-scalar and iso-vector mesons 



(79) 



We investigate this relation with our numerical results 
later. At first, we need to solve the eigenvalue equation 
(61) numerically to obtain wavefunction of the nucleon 
fh.R and the nucleon mass tub- It is useful to define the 
dimensionless variables w = Mkkw, U = U/Ukk, and 
z = z/Ukk which are related together as 



w 



dz 



o [i + z 2 p 2 ./, y r r - - i 



dU\ 



U 



(80) 



We can rewrite the m b (w) in terms of the new dimen- 
sionless variables as 



m b (w) ~ mf ] ■ U = M KK m b (w), 



(81) 



where 



. . v /372 4 7 r 2 /i4 (2ira') 2 R N 

m b (w) = y~L '- N c U(w). (82) 
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After rewriting the equation (61) in terms of w, we obtain 



[-d% - dwm b {w) + (m b (w)) 2 ] f L (w) 



m B 



M. 



KK , 



(83) 
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The key idea for using dimensionless variable is that the 
functions /z,(u>) does not depend on the scales further. 
Now, we use the shooting method again to solve the 
above equation numerically and find the Jl{w) and its 
eigenvalue tub /Mkk- In order to do the numerical com- 
putation, we assume N c = 3 for realistic QCD. Also as 
it is mentioned in the previous section, we fix the value 
of Mkk = 0.395 GeV to have the pion decay constant 
f v = 0.093 GeV. We obtain the various coupling by cal- 
culating the integrals (66) and (71) and compare some of 
our results with the results of SS model [23] in Table 2. 



TABLE II: . Numerical results for Axial and vector meson 
couplings in the non-critical holographic model of QCD. The 
value of vector couplings are compared with the SS model 
results[23]. Superscript (a) and (b) are stand for our results 
and the SS results, respectively. 



k 


(k) 

9A,min 


(k) 
9 A ,ma 9 


(k).a 
9v,Ln 


<k),b 
9 V ,miu 


<k).a 
9 V ,ma, 


(k),b 
9v,maq 





1.16 


1.86 


8.30 


5.933 


-1.988 


-0.816 


1 


1.07 


1.44 


1.6488 


3.224 


-6.83 


-1.988 


2 


0.96 


0.862 


1.9 


1.261 


-7.44 


-1.932 


3 


0.67 


0.14 


0.688 


0.311 


-4.60 


-0.969 



Also, we obtain the axial couplings as 

9A,mag — 1-582 , gA,min — 



(84) 



while in the previous analysis done [23] using the SS 
model, it is obtained 



9A,r. 



9 A,; 



0.13 



(85) 
0.7 and 



which if we consider N c = 3, then gA.mag 
9 A = 0.83. It should be noted that the higher order 1/N C 
corrections can be used to improve this result. The lattice 
calculations indicated that higher order 1/N C corrections 
are suppressed. Our result are a good approximation of 
the experimental data at leading order g^ p — 1.2670 ± 
0.0035. 



V. NUCLEON-MESON COUPLINGS 

In the previous section, we demonstrated that the 
nucleon-nucleon interaction can be interpreted as me- 
son exchange potential and obtained the coupling con- 
stants in different channels using a non-critical holo- 
graphic model of QCD. Now we aim to calculate the iso- 
spin dependent and iso-spin independent pseudo scalar 
meson, vector meson and axial vector (pseudo-vector) 
meson couplings to the nucleon. The couplings of vec- 
tor meson to the nucleon can be related to the vector 
couplings as 

AT (k), singlet at (*0 

jv c 9y _ Jv c yy,m in 

9uj( k )NN — 



2 2 

(k), triplet ( k ) , ( fe ) 
gy _ 9v,min ' 9v,mag 

9 p Wnn — 2 — — ' 2 



(86) 



For the axial vector meson such as which is an iso- 
vector meson, there is 

(k), triplet ( k ) , ( k ) 

9A 9 A.min ' 9 A,mag ron\ 

9aWNN = ^ = 2 ' ' 

and for the iso-scalar axial vector meson we have 



at Ak), singlet jy O) 



9f( k )NN = 



The pseudo-scalar mesons such as pion it and 77' couple 
to the nucleon as 

9ir( k )NN TiT _ 9 A n/r 9 Amin + 9 A.mag 

mkk = —^— r — m K k = — M KK , 



2m,N 

g,f( k )NN 

2m n 



M 



KK 



2/tt 

N c g s 2 nglet 
2/. 



M. 



KK 



2 h 

N c 9A,m 

2U 



Mkk- 



(89) 



We calculated all of the leading order meson- nucleon cou- 
plings numerically and show them in Table. 3. We 
compare our results with the modern phenomenological 
nucleon-nucleon interaction models predictions such as 
the AV 18 [8], the CDBonn [7], the Nijmegen(93) [6] 
and the Paris [5] potentials.lt is necessary to mention 
here that the components of these models are very differ- 
ent in strength, and if parameterized in terms of single 
meson exchange give rise to effective meson-nucleon cou- 
pling strengths, which also are similar. In Table 3, we use 
the value of effective meson-nucleon couplings for these 
potentials. 



TABLE III: . The values of different effective meson-nucleon 
couplings in the phenomenological interaction models [54], SS 
model [23] and our model. 



g 


V18 


CD - Bonn 


Nijmegen (93) 


Paris 


SS model 


Our model 


9a« 


9.0 


9.0 


9.0 


10.4 






9a 


9.0 


11.2 


9.8 


7.6 






9* 


13.4 


13.0 


12.7 


13.2 


16.48 


15.7 


9v 


8.7 


0.0 


1.8 


11.7 


16.13 


0.0 


gu, 


12.2 


13.5 


11.7 


12.7 


12.6 


12.45 


9p 




3.19 


2.97 




3.6 


3.15 


9a> 










3.94 


1.51 


9p 












1.74 



The isospin independent scalar potential component 
in the interaction models is due to a single scalar meson 
exchange interaction. The mass of lowest scalar meson is 
not established well [55], but in the phenomenological 
NN interaction models it is typically taken to be of order 
500 — 700 MeV. By considering ms to be roughly 600 
MeV, the effective scalar meson coupling constant for 
these interaction models differ from 7.6 (Paris) to 11.2 
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(CD-Bonn). It should be noted again that these values 
are the effective couplings. In fact two of these models 
don't contain any scalar meson in their parameterized 
forms. In our holographic model based on the non-critical 
string theory, there is no scalar interaction term too. The 
isospin dependent component of scalar potential can be 
interpreted as a scalar spin 1 meson exchange. In N c 
counting, these component are of order 1/N C , so they are 
weaker than the isospin independent scalar mesons by of 
order of their magnitude. The effective values for the 
lowest scalar meson (a (980)), range from 9.0 to 10.4 in 
the various phenomenological NN interaction models. 

The vector component of phenomenological NN inter- 
action models considered here, is equal to the scalar com- 
ponent with the minus sign which means the strength of 
a (980) exchange interaction is equal to the exchange of 
p(770) in the large N c limit. Indeed the vector interaction 
arises from a strongly ui exchange interaction for isospin 
independent component and a weaker p meson exchange 
for isospin dependent component. The value of effective 
uj— nucleon coupling range from 11.7 (Nijmegen(93)) to 
13.5 (CD-Bonn), while in the original version of these 
models this value varies from 10.35 (Nijmegen) to 15.85 
(Bonn) [56]. In the analysed done in the SS model, by 
considering M KK = 940 MeV, N c = 3 and A = 17 which 
imply rriN — 1.93 Mkk, 9oj evaluated 12.6. We also ob- 
tain the value of g u = 12.45 which is in the range of the 
values anticipated from the phenomenological potential 
models. We have used the iV c = 3, also Mkk = 395 
MeV. The obtained value for the nucleon mass in our 
model is roughly 920 MeV which is very realistic and close 
to the familiar nucleon mass. The isospin dependent com- 
ponent of vector potential which arises from a p meson 
exchange potential is roughly three times weaker than 
the isospin independent component. In a chiral quark 
model, it is expected to have g u = 3g p , but the value 
of the 1Z = .9^/3 g p differs from one in the above phe- 
nomenological interaction models. It is 1.66 for the CD- 
Bonn, 1.5 for the Nijmegen and 0.77 in the Paris model. 
This ratio is about 1.2 in the SS model and equals to 
1Z = 1.33 in our model. Actually, the NN phase shifts 
uniformly require a larger 1Z than the chiral quark model 
prediction which is a mystery. But in the potentials come 
from the holographic QCD model, it can be explained by 
the contribution of the magnetic coupling in the vector 
channel. 

The Nijmegen(93) and CD-Bonn models don't contain 
any single axial vector meson exchange, so there is no 
axial vector interaction in their structures. The phe- 
nomenological AV 18 and Paris potentials, predict a small 
value for the axial vector interaction. In our calculation, 
the <X(i)— nucleon coupling equals to 1.51 and the value 
of /(i)— nucleon coupling is obtained 1.74. 

The isospin independent pseudoscalar interaction 
comes from a 77 meson exchange. This component isn't 
well constraint by NN scattering data and the phe- 
nomenological interaction models give extremely differ- 
ent values for this component ranged from 1.8 (Ni- 



jmegen(93)) to 11.7 (Paris). While analysis of other ob- 
servable such as rj— meson photoproduction suggest that 
the coupling constant value shouldn't go beyond the 2.2 
[57]. In our calculations, g v is in the leading order. 
Despite the isospin independent component, the isospin 
dependent component is strong and spread on a long 
range, all of the models considered here have a main 
component for this interaction which is a pion exchange 
potential. The values of the pseudoscalar pion-nucleon 
coupling constant g^ vary from the 12.7 to 13.4 effec- 
tively. In our calculation, g^ is evaluated 15.7 while in the 
calculations based on the SS model was obtained 16.48 
[23]. 



VI. CONCLUSION 

In this paper we have used the non-critical holographic 
duality of gauge theories in the background of the near- 
extremal AdSe. We obtained the mass scale of the model 
by comparison of the pion field action with the usual 
Skyrme model action. We show that the size of baryon 
is of order one in respect to the t' Hooft coupling but is 
still smaller than the scale of the dual QCD, so we in- 
troduced an isospin 1/2 Dirac field for the 5 dimensional 
baryon and wrote an effective action for it. Then, we 
reduced the 5D action to 4 dimension by using a mode 
expansion for the baryon field along with the mode ex- 
pansion for the gauge field on probe branes and attained 
the nucleon-nucleon interaction in terms of meson ex- 
change potentials. 

In our analysis such as the critical ones, the NN po- 
tential involves only the pseudoscalar, vector and axial 
vector meson interactions. In fact these meson exchange 
potentials play a special role in producing the nucleon- 
nucleon potential. The long-range part of the nucleon- 
nucleon potential (r > 3/to) is mostly due to the one 
pion exchange mechanism which is the strongest compo- 
nent among the iso-spin dependent components. In the 
phenomenological potential models, the isospin depen- 
dent pseudo-scalar meson exchange potential is of order 
N% while the isospin independent component is of order 
N c . So it is expected that the 77- meson exchange poten- 
tial would be very weaker than the the 7r-meson exchange 
potential. We obtained the g^ = 15.7 and g v = 0.0 in 
the leading order that is consistent with the phenomeno- 
logical interaction models. 

Isospin independent scalar mesons are responsible for 
the attractive interaction in the intermediate range of 
the potential (0.7 < r < 2fm). This components are the 
main reason for the nuclear binding. Also in the phe- 
nomenological interaction models, the strength of this 
interaction is equal to the vector meson exchange with 
a minus sign. In fact the radial shapes differ consider- 
ably at short distances, ranging from attractive to re- 
pulsive. Some of this potential models just involve the 
scalar meson exchange and another contain the vector 
meson exchange interaction. We consider the vector me- 
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son exchange in our analysis which produce the strong 
short-range repulsion. Exchanging the vector meson p 
can explain the small attractive behavior of the odd- 
triplet state. 

We compared our results with the available values of 
coupling constant predicted in four modern phcnomcno- 
logical interaction models (Nijmegen(93), Paris, CD- 
Bonn and AV 18 models) in Table. 3. The remarkable 
point here is that we calculated all of the meson-nucleon 
couplings whereas in the phenomcnological interaction 
models these values fit to the NN scattering data and as 
it is indicated from Table 3. the value of coupling con- 
stants are widely different in various interaction models. 

We believe that our noncritical holography model is 
more reliable than the critical SS model to study the 
nucleon-nucleon interactions for some reasons. At first, 
the size of baryon was proportional to 1 / \/XN c in the 
SS model and becomes zero at large 't Hooft coupling. 
While in our calculation we show that the size of baryon 
is of order one. So, our model doesn't suffer from the 



zero size of baryon. 

Also, the nucleon mass was about 1.93 M KK in the SS 
model. The mass scale was obtained Mkk = 940 MeV 
to describe the meson spectrum. So the nucleon mass 
was about 1.8 GeV there. It caused some inconsistency 
in analysing baryon in the SS model, because the scale 
of the system associated with the baryonic structure is 
roughly half the one needed to fit to the mesonic data 
[36] . But in our analysis the mass of nucleon is obtained 
roughly 920 MeV which is very close to the neutron (pro- 
ton) mass. In addition our noncritical calculations can 
describe the meson-nucleon couplings successfully at least 
at the leading order. 

We can consider the derivative couplings to the 
nucleon-nucleon interaction and obtain the contribution 
of other mesons to the potential. Also, the couplings can 
be computed with more accuracy. Moreover, we can use 
the obtained potential here to study the nuclear proper- 
ties such as the NN scattering data or nuclear binding 
energies. We leave them here as the future works. 
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